TOWARDS THE KPP-PROBLEM AND logt-FRONT 
SHIFT FOR HIGHER-ORDER NONLINEAR PDES III. 
DISPERSION AND HYPERBOLIC EQUATIONS 



V.A. GALAKTIONOV 

Abstract. Some aspects of extensions of ideas of Kolmogorov, Petrovskii, and Piskunov 
(1937) |13) on travelling wave propagation in the reaction-diffusion equation 

Ut = + u{l — u) in M X IR+, Uq{x) — H{—x) = {1 for x < 0; for a; > 0}, 

are discussed. The present paper continues the study began in [21 [TU] for higher-order 
parabolic semilinear and quasilinear bi-harmonic equations such as 

Ut ^ -u^xxx + uil - u), ut^ ~{\u\'^u)xxxx + u{l~u) (n>0), etc. 

Here, higher-order dispersion equations such as (Dx = ^ and Df = 

Ut = —D]^u + 7i(l — u) and up to D^u — —D^^u + u{l — u), etc. 

are studied. Some features of KPP-like results are also shown to exist for semilinear 
dispersion-parabolic equations such as 

Uttt = ~D}^^u + u(\ — Uttttt — D}^u + u(l — u), and others, 

and for pure hyperbolic ones 

utt = -Uxxxx + u{\ ~ u) and up to utttt = ~ D]^ u + u{l ~ u) , etc. 

As an example, we also treat a quasilinear ViyEj ut — — -D^^(|u|"u) -f w(l — u), with n > 0. 
Two main questions are: (i) existence of travelling waves via any analytical/numerical 
methods, and (ii) their stability and derivation of the log ^-shifting of moving fronts. 



1. Introduction: the classic KPP-problem and other higher-order PDE 

MODELS 

1.1. The classic KPP-problem of 1937. The classic KPP-problem [I3] (1937) 
(1.1) Ut = Uxx + u{l — u) in M X M+, u{x,{)) = uq{x) in M, 
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with the step (Heaviside) initial function 
(1.2) uo{x) = H{~x) = 



1, X < 0; 
0, X > 0, 



consists of studying large-time convergence of the solution of the Cauchy problem ( II. ip . 
(11. 2p to the unique minimal travelling wave solution, with the minimal speed Aq = 2 of 
the standard form 



u^{x,t) = f{y), y = X - Xot, where 
-Aof = /" + /(I - /), 2/ G M; /(-oo) = 1, /(+oo) = 0. 



(1.3) 

The KPP paper [13] contains a number of pioneering remarkable results, which founded 
several new directions of modern nonlinear PDE theory. For further use, it suffices for us 
to refer to a survey and key references in [9]. 

Let us state the main result of |T3]. The convergence to the minimal TW (I1.3P was 
performed in the TW moving frame, proving that the TW front moves like 

(1.4) Xf{t) =2t- g{t) as t ^- oo, with g{t) = o{t), 
where the front location Xf{t) is uniquely determined from the equation 

(1.5) u{xf{t),t) = \ for all t > 0. 
Then the convergence result of [13] takes the form: 

(1.6) u{xf{t) + y,t) f{y) as t +oo uniformly in y G M. 

In 1983, Bramson [1], using probabilistic techniques, proved that there exists unbounded 
log t-shift of the moving TW front 

(1.7) ^(t) = A;logt(l + o(l)), withA; = |, 

Therefore, (11. 7p implies eventual, as t — )■ +oo, infinite retarding of the solution u{x,t) 
from the corresponding minimal TW, thought the convergence (ll.6p takes place in the 
TW frame. 

1.2. KPP-like problem to higher-order semilinear and quasilinear parabolic 

PDEs. We dealt with such semilinear reaction-diffusion PDEs in [9J, where the main 
basic model was the semilinear bi-harmonic equation, i.e., a fourth-order semilinear heat 
equation (SHE-4) 

(1.8) ut = -u^ccxx + u{l -u) in M X M+. 

The corresponding TW with the speed of propagation A is then governed by the following 
fourth-order ODE: 

(1.9) u4x,t) = f{y), y = x-Xt =^ _A/' = -/^ + /(I - /), 
with the singular boundary conditions at infinity: 

(1.10) f{y) — and f{y) — )■ 1 as ?/ — )■ ±oo "maximally" exponentially fast. 
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In particular, we have found a "maximal speed Amax = 1.27148... such that 

(1.11) TW profiles f{y; A) exist for all < A < Amax, and nonexistent for A > A^ax- 
In ^Q\, we extend some of the above results to the quasilinear KPP-4n problem for 

(1.12) ut = -{\u\'u)^^^^ + u{l-u) in RxM+, 

where n > is a parameter, as well as to some other parabolic equations. 

In a similar manner, we studied in [9] the semilinear tri-harmonic equation (SHE-6): 

(1.13) ut = u^xxxxx + u{l -u) in M X ]R+, u{x, 0) = H{-x) in M. 

The corresponding TW with the speed of propagation A is then governed by the following 
fourth-order ODE: 

(1.14) u.{x,t) = f{y), y = x-Xt =^ -A/' = /(^^ + /(I - /), with ([HQ]). 

The (II. lip remains true, with Amax = 2.12110... . 

In general, in [5], we presented some numerical evidence on existence of various TWs 
and there properties for semilinear parabolic 2mth-order PDEs (SHE-2m) such as (here, 
D = ^) 

^ dx' 

(1.15) ut = {-lT+^Df'u + u{l-u) [Dx = j^), with the ODEs 

(1.16) u,{x,t) = f{y), y = x-Xt =^ _A/' = (-1)^+1/(2-) + /(i _ /) 

(plus (11.101) ). and rather sharply estimated Amax = Amax("^) > for m = 3,4,5, i.e., up 
to the tenth-order parabolic equation as in fll.l5p . 

1.3. Results I: dispersion PDEs (Section [2|) . In the present paper, we will deal with 
higher-order dispersioi^ equations 



1.17) 


Ut = 


-Dl\ + u{l 


-u) = 


-xf' = 


_/(n) + j(l_ 


/); 


1.18) 


Urn = 


-Dl}u + u{l 


-u) = 


^ -x'r- 


= _/(!!) + /(I. 


-/); 


1.19) 


Uttttt = 


-Dl\ + u{l 


-u) = 


-x'f'^ 




-/) 


1.20) 


Dju = 


-Dl\ + u{l 


-u) = 






-/) 



^Here, we use a PDE classification, associated with some a priori bounds admitted by the principal 
linear differential operators by multiplying by appropriate time derivatives D^u in the L^-metric. There- 
fore, while, for most of parabolic and hyperbolic equations, such a classification coincides with the classic 
rigorous Petrovskii's one P2] (recall: "parabolic in Petrovskii's sense"), for others, our classification can 
be different and uses a "more applied" understanding of dispersion and related phenomena. Anyway, we 
do not think that, for truly higher-order (in both x and t up to 11th or 12th orders) semilinear PDEs, any 
formal classification may somehow essentially help to understand the nature of TW patterns obtained 
below (note that such "evolution" patterns have been obtained even for obviously elliptic PDEs, for which 
a standard evolution interpretation makes no sense, due to Hadamard's example of an ill-posed Cauchy 
problem) . 
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(1.21) Dlu = -D':u + u{l-u) =^ _a9/(9) = -/(") + /(I-/); 

plus singular boundary conditions f ll.lOp . We have chosen rather higher-order dispersion 
operator D]^ on the right-hand side in order to avoid a kind of "temptation" to rely 
on linearized analysis, which was heavily used in [Qj § 2] in studying the lower-order 
bi-harmonic KPP-problem for fll.Sp . 

1.4. Results II: dispersion-hyperbolic PDEs (Section [51). We next consider exis- 
tence of TW solutions of higher-order dispersion-hyperbolic equations 



(1.22) 


Utt = 


-Dl'u + u{l - 


u) = 


> AY' = - 


-/^'^ + /(!-/); 


(1.23) 


Utttt = 


-Dl\ + u{l - 


u) = 


> aV^'^ = 


-/^'^ + /(!-/); 


(1.24) 


Utttttt = 


-Dl^u + u{l- 


-u) = 




-/("^ + /(!-/); 


(1.25) 


D'tu = 


-Dl\ + u{l - 


u) = 


> aV^«) = 


-/^'^ + /(!-/); 


(1.26) 


D]\ = 


-Dl}u + u{l - 


u) =^ 


^ ^10 j(10) ^ 





plus singular boundary conditions fll.lOp . 

1.5. Results III: dispersion-parabolic PDEs (Section [4]). We next consider, in the 
KPP setting, the following equations, which, in view of certain estimates, can be con- 
sidered as dispersion-parabolic (the signs in front of avoids "backward parabolic" 
features) models 



1.27) 


Uttt = 


-Dl\ + u{l 


-u) = 




_/(io) + J(l _ 


-/); 


1.28) 


Uttttt 


= Dl\ + u{l 


-u) = 




/(lO) + J(l _ 


/); 


1.29) 


Dlu = 


-Dfu + u{l 


-u) = 




_/(10)+J(l 


-/) 


1.30) 


D'tu 


= Dfu + u{l 


-u) = 


=> -aV(9) = 


/(lo) + ;(i _ 





with the conditions fll.lOp . 

1.6. Results IV: higher-order hyperbolic PDEs (Section E]). Finally, we consider 
four purely hyperbolic higher-order equations: 



1.31) 


Utt = Dl\ + u{l 


-u) = 


> x'f" = 


= /(^°) + /(l-/); 


1.32) 


Utttt = -Dl% + u{l 


~u) = 


> aV^'^ 


= -/(i°) + /(l-/): 


1.33) 


D^u = Dl% + u{l 


-u) = 


> aV^'^ 


= /('°) + /(!-/); 


1.34) 


D^u = -Dfu + m(1 


-u) = 


> aV^'^ 


= + /(I - /) 
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with the conditions fll.lOp . We also present some TW patterns for the corresponding 
elliptic PDEs such as 

(1.35) uu = -Dl\ + u{l-u) =^ aV" = -/^'°^ + /(!-/); 

(1.36) Uuu = Dl% + u{l-u) =^ x'fi^) = fm + f(i-f). 

(1.37) Dlu = -D^^u + u{l-u) =^ x'f^'^ = -f^''^ + f{l-fy, 

Clearly, as evolution equations, these lead to unstable (and ill-posed, in Hadamard's sense, 
1906) problems, but as certain special TW patterns, such solutions make sense, though 
are highly oscillatory, as one can expect. 

1.7. Results V: KPP-(10,11) and KPP-(11,12) (Section[6]). These are most exotic 
KPP-models under consideration. As an example, we consider two of such PDEs, where 
the order in t exceeds the order in x: 

(1.38) Dl^u = -Dl% + u{l-u) =^ -A"/(ii) = + /(1-/), 

(1.39) Dru = -Dl\ + u{l-u) =^ AiV^i2) = -/(ii) + /(l-/), 
with conditions fll.lOp . 

1.8. Results VI: semilinear eleventh-order in time PDEs (Section [7]). We end 

up with the following equations, which were not treated above and are induced by fll.38p : 

(1.40) Dl\ = -u, + u{l-u) =^ -A^V^"^ = -/' + /(!-/), 

(1.41) Dl'u = -u,, + uil-u) =^ _AiV(ii) = -/" + /(1-/), 

(1.42) Dl\ = -u,,, + uil-u) ^ -Aii/dD = -/'" + /(I-/), 

(1.43) Dl'u = -u,,,, + ^.(1 -u) =^ -A^V^ii) = -/(^) + /(I - /). 

Note that, according to a priori bounds, the models f ll.4ip and f ll.43p belong to the 
parabolic type, while fll.401) and fll.43p can be treated as dispersion ones. 

1.9. Results VII: a quasilinear dispersion KPP— (11,1) problem (Section [8]). As 
in [in] for parabolic KPP-problems, we claim that many present results admit extensions 
to quasilinear PDEs. To this end, we treat the following quasilinear KPP-(11,1) (cf. 
f ll.l7p ) problem 

(1.44) u, = -DlWuru) + uil -u) =^ -\f = -d/r/)^") + /(I - /), 

where n > is a fixed parameter. The main distinguished feature of degenerate equations 
as in fll.44p is that the TW profiles f{y) have finite interface at some finite yo > 0, so we 
need first to explain the local (periodic) structure of such solutions as y ^ . 
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1.10. The origin of logt-front shift and general goals of the paper. In Section [9l 
using a general third-order in t semilinear model, we show, via a kind of an "affine centre 
subspace expansion" , how log t- front shift can appear for t ^ 1 in convergence to a TW 
pattern. 

Concerning our general classification, using the KPP-setting, we refer to such problems 
as to the KPP-(A;, /), where k stands for the order of the differential operator in x and / 
for the order of the derivative in t. 

Thus, for several KPP-(A;, /) problems, with k > 3 and / > 1, the main questions to 
study, here, as in P [10], are: 

(I) The problem of TW existence: existence of travelling waves by using analyt- 
ical/numerical methods, 

(II) The problem of "maximal" speed: defining, in a natural sense, the a;-limit 
set uj{H) of the properly shifter orbit (11.61) . to discuss whether or not, at least, in some 
particular KPP-problems, 

(1.45) uj{H) = {/(■; A) : A G A}, where A C R is bounded. 
Two main questions arise: 

(1.46) what is the maximal speed Amax = supA, and whether A = {Amax}? 

Those questions are "remnants" of the KPP setting. It turns out that, unlike in [13], for 
several parabolic [9l [10] and dispersion models 

(1.47) A = (0, Amax), so that A^ax = sup A, Amax ^ A, and A^in = ^ A. 

However, for some PDEs under consideration, there exist stationary patterns correspond- 
ing to A = 0. Therefore, for such a variety of semilinear PDEs of different orders, any 
clear general conclusions on the above speed sets are illusive, and each problem might 
have some individual features. 

(III) The logt-SHlFT problem: studying the stability of the TW /(y; Aq) and deriva- 
tion of the log t-shifting of the moving front in the problem (II. ip . (11.21) . connected with a 
kind of an "affine centre subspace behaviour" for the rescaled equation (Section [9]). 

Overall, as in P [ID], we expect that the log t- shifting phenomenon is quite a generic 
property of many nonlinear KPP-type problems regardless their particular types. 

Note that, here, we are not able to solve or even discuss one of the main problem about 
the TW velocity Aq (or velocities? - the oj-limit set of the rescaled orbit might include 
a connected continuous curve of profiles {f{-,X), A G A}), which appear in the PDE 
setting with the Heaviside initial data (II. 2p (cf. (II) above). Numerically, this would 
require a full set of hard PDE numerical experiments, which we do not perform in such 
a generality. The author believes that using such a full-scale of both the ODE and PDE 
numerics is too exhaustive and inevitable moves the research into a pure numerical area, 
where some important mathematical aspects of the KPP-ideology, which comprise the 
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main goal of the present study, could be lost would be essentially reduced to numerical 
aspects to appear for sur^. 

2. Dispersion equations of various orders in t 

As in f9] for PDEs (11 .Sp . fll.lSp . and f ll.lSp {m = 4,5), for new classes of dispersion 
PDEs fll.17p - fll.2ip . it is convenient to present numerical results, which directly show the 
global structure of such TW profiles to be, at least partially, justified analytically. A 
more detailed description of using the bv4p4c solver of the MatLab is given in [9l § 2]. 
However, it is important to note that, as the initial data for further iterations, we always 
took the Heaviside function as in fll.2p . This once more had to help us to converge to a 
proper "minimal" profile (indeed, there are many other TW profiles), though, of course 
this was not guaranteed a priori. We keep this rule for all other KPP-{k, I) problems of 
interest, including the hyperbolic ones fll.3ip - fll.34p and more higher-order ones, where 
initial velocity, acceleration, etc. were taken zero. 

Thus, numerically, we observe existence of TWs for all the PDEs fll.l7l) - fll.21l) . In 
Figure [H we show the TW profiles for the first-order (in t) dispersion equation (11.171) for 
the speeds A = —1, 0.5, and 1. In addition, a stationary solution (denoted by the dash 
line in Figured]) exists for A = 0, hence, satisfying the ODE (to be referred to later on a 
few times) 

(2.1) /(!!) = /(I-/) in M, plus dnoD. 

Existence (but not uniqueness) of a solution to (12.11) (and of ones of arbitrary odd order 
j{2m-i) ^ with any m > 2) is well known for a long time; see [21 [T3j. Uniqueness of such 
stationary solutions of equations like (12. ip of various odd orders was observed numerically; 
see [7j, where further related references were available. For convenience and for future 
use, in Figure [31 we present the solution of this stationary (A = 0) problem (12.11) . 

Figure [2] shows a quite oscillatory TW profile for A = —3. 

Further numerical experiments confirm the following A-range of existence of such TW 
profiles for ffLTTD (cf. ffLTTj) ): 

(2.2) KPP-(11,1): 1.2< An,ax< 1.3. 

The TW profile for the existence parameter A = 1.2 is shown in Figure [H Numerical 
integration is performed on a smaller interval [—60,60], that, in view of small oscillatory 
behaviour of f{y) in Figure [H, cannot cause any problem (though the length of the interval 
might affect this critical speed value). Moreover, it is seen that, as A — )■ A"^^^^, the profile 
f{y) tends to be non-oscillatory at all, so that this Amax, characterizes the case when the 
characteristic equation of the linearized operator admits complex roots with vanishing 

^In other words, the author's position here as follows: (ODEs) numerics, indeed, helps to understand 
various features/aspects of modern nonlinear higher-order PDE theory, but should be used carefully and 
balanced, i.e., not shadowing (and not replacing!) any mathematical ideas/methods/tools/ideology/etc. , 
which could appear from and connected with treating reliable numerical results (though the above is 
more relevant to [9] than to the present paper). 
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f(y) KPP-(11,1): -;>.f'= - f'"' + f(1-f):TW profiles forX = -1,0, 0.5 

I.4F 

X=-1 




Figure 1. A TW profile f{y) satisfying tfie ODE in ffLTTj) for A 



f ^ KPP-(1 1 , 1 ) : -X f = - f<" ' + f(1 -f), TW profile for X = -3 




-200 -150 -100 -50 



50 100 150 200 



Figure 2. An oscillatory TW profile f{y) satisfying the ODE in (11 .17^ 
and (fTTOj) for A = -3. 



imaginary part. Then this changes the dimension of the asymptotic bundles and makes 
matching impossible; cf. Theorem 2.1 in [Qj § 2]. 

For equations (I1.18p - (ll.2ip . TW profiles, with various negative and positive A's, are 
presented in Figure |5l In all the cases, there exists also the stationary profile for A = 
satisfying (12. ip . 




-50 50 

y 

Figure 3. The stationary TW profile f{y) satisfying the ODE in f ll.lTp 
for A = 0. 




y 

Figure 4. A TW profile f{y) satisfying the ODE in ffTTTll for A = 1.2 ^ A^ 



3. Dispersion-hyperbolic PDEs 

For equations (11.22p - fll.25p . TW profiles, with various speeds A's, are presented in 
Figure [61 In the last case (c), i.e., in the KPP-(11,8), numerics reveal existence of Amax 
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KPP-(1 1 ,3): -X^ r = - f + f(1 -f): TW profile fori = -0.6 and 0.6 




KPP-(1 1 ,5): -r" f'^'= - f<"' + f(1 -f): TW profile forX = -0.7 and 0.7 



(a) (fLT8)) 

KPP-(1 1 ,7): -X'r = - + f(1 -f): TW profile for>. = -0.8 and 0.6 





(b) (HHl) 

KPP-{1 1 ,9): -7? f'''= - f'" ' + f(1 -f): TW profiles fori = -0.9 and 0.9 




-40 -30 -20 -10 10 20 30 40 



(c) i^rm 



(d) dm]) 



Figure 5. TW profiles for dispersion equations (fTTT8|) - pr2T]) . 



satisfying 

(3.1) 1.0443 < A„,ax < 1.0445. 



When approaching X^^^^, the TW profiles f{y) remain essentially oscillatory for y 1, 
so nonexistence for slightly A > Amax is not related to changing of the dimension of the 
linearized bundle (cf. |i9] § 2.5] for the KPP-(4,1)), but means the impossibility of the 
corresponding "nonlinear matching" of such bundles. 

For the most exotic KPP-(11,10) equation fll.26p . we present TW profiles for A = 
(the stationary equilibrium), A = 1, and A = 1.1 in Figure [71 Even for A slightly larger 
than 1, the profiles get very oscillatory about the equilibrium 1 as 2/ — >■ — oo; cf. A = 1.1 
in Figure [71 
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f(y) 



KPP-(11,2): r f"=-("" + f(1-f):TW profiles for). = 0,0.5,1 





[ 


>.=0.5 


■ ^^^^ ^ 
J' 








-50 -40 -30 -20 -10 10 20 30 40 

y 



(a) (HHa), A = 0,1,0.5 

f{y) KPP-(11,6): x'' f<'"=-f''" + f{1-f):TW profiles fori = 0.6 and 1 




f(y) 



KPP-(1 1,4): >.''f''"= - f'"' + f(1-f): TW profiles fori = 0,0.6,1 



-^Jf 




1=0.6 ^^"^ 









-40 -30 -20 -10 10 20 30 40 

y 



(b) (fTMl) . A = 0,0.6,1 

f(y) I^PP-I^ 1.8): ).''f''''= - f'"' + f(1-f): TW profiles for). = 0, 0.6,1 ,1 .04 



— -4 


; _ >.=o, 0.6 


1 ».=i k 


^^^^^^^ 







-40 -20 



20 40 60 

y 



(c) (1123, A = 0.6,1 (d) dOSl), A = 0.06, 1, 1.04 

Figure 6. TW profiles for dispersion- hyperbolic equations (|1.22p - (|1.25p . 



4. Further dispersion-parabolic equations 



For equations ( 11.27p - p.30p . TW profiles, with various speeds A's, are presented in 
Figure [Hi For equations fll.281) and fll.301) . for A = 0, the stationary equation 



(4.i: 



/(10)=_J(1_J) in 



admits a solution with a periodic behaviour as y — > +oo; see Figure [9l 

Similarly, by symmetry, the stationary equations f ll.27p and (11. 290 admit analogous 
stationary profiles that are oscillatory about 1 as y — ?■ — oo. 



5. Higher-order hyperbolic equations and elliptic patterns 



5.1. Hyperbolic equations. For equations fll.3ip - (ll.34p . TW profiles, with various 
speeds A's, are presented in Figure [101 
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-50 



y 

Figure 7. TW profiles for tlie KPP-(11,10) problem ([L26]). 



Note that, for such hyperbolic problems, the behaviour as y — )■ +oo often becomes 
non-decaying oscillatory (as a feature of hyperbolic flows), without a decay or with a very 
slow decaying algebraic envelope. However, as seen from some figures above, several TW 
profiles decay at +oo, and hence satisfy the standard KPP setting. 

Thus, here, we fix operators that are tenth-order in x to avoid any questions on a 
possibility of a reliable and rigorous local and/or global ODE analysis similar to that 
performed in [9, § 2] for the parabolic KPP-4 problem. However, the asymptotic study 
of these ODEs as y — )■ +oo can be performed justifying that the dimensions of stable 
bundles as y — > ±oo well-correspond to the existence of TW profiles in the sense of 
a multi-parameter shooting. Nevertheless, any rigorous proof of existence of such f{y) 
remains hopeless and represents an open problem, as the existence of a heteroclinic path 
between the equilibria and {1,0,..., 0} in the tenth-dimensional phase space occurring 
for ODEs in (fCTjl -f lLM]) . 

Any stability analysis of such TWs in both hyperbolic and dispersion cases leads to very 
difficult spectral problems for pencils of non self-adjoint operators, which we are not going 
to treat here; see further "spectral" references and examples for hyperbolic/dispersion 
equations in [SI § 8.4]. 



5.2. Elliptic patterns. Though the corresponding elliptic equations (ll.35p - (ll.37p do not 
admit any evolution setting, in Figure [TTl we present elliptic patterns, which are always 
highly oscillatory and "almost periodic" either as y ^ — oo or y — )■ +oo. 
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f(y) 



KPP-(10,3): -).f=-f'"' + f(1-f):TW profile for). = -0.5 




r 



KPP-(1 0,5): - ).^l'^'= f''°'+t(1-f), TW profile for ^ = 1 



-250 -200 -150 -100 




-30 -20 -10 10 20 30 40 50 



(a) nrm . x ^ -0.5 

f(y) KPP-(10,5): -X^f'^'s f'''"+f(1-f), TW profile for). = 0.5 




f(y) 



(c) (fLMI . A = 0.5 

KPP-(10,7): -).'f'"=- f"°' + f(1-f), TW profile for>. = -1 




-50 -40 -30 -20 -10 10 20 30 40 

y 



(b) (HHl), A = 1 

KPP-(1 0,5): - f °'+f(1 -f), TW profile for X = 1 ,2,3,4 




(d) (fLMI . A = 1,2,3,4 

f(y) KPP-(10,9): -I'f^L f'"'+f(1-f), TWprofile for>.= 1 





-30 -20 -10 10 20 30 40 50 60 

y 



(e) (11211, A = -1 (f) (ESni), A-1 

Figure 8. TW profiles for dispersion-parabolic equations (|1.27p - (|1.30p . 
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KPP-(10,5): f ' = f(1-f): stationary TWprofile for?. = 



71- 




1 


m 




_ 










-20 -10 10 20 30 40 

y 


50 60 



Figure 9. A stationary, A = 0, solution of fl4.ll) with a periodic behaviour 
at +00. 

6. KPP-(10,11) AND KPP-(11,12) 

In Figure [121 we present TW profiles for the KPP-(10,11) problem fll.381) for A = 0.2 
(a) (after a proper reflection, this profile is very close to the stationary one in Figure [H]) 
and for A = -1, -1.1 (b). 

In Figure [131 we show TW profiles for the PDF (11.391) for A = 1 and 0.8. 



7. Some semilinear PDEs that are eleventh-order in t 

For equations (ll.40p - (ll.43p . which all of different types, the TW profiles look rather 
similar and are presented in Figure [T^ Note that, in all the cases, we fix negative A = — 1. 
Choosing positive speeds A = +1 always led to highly oscillatory behaviour at the right- 
hand end point of the interval of integration and no convergence was observed. 

Note also that the first equation (11.401) in this list admits obvious (and well known) 
explicit stationary solutions for A = 0: 

(7.1) (11301), A = 0: f = -/(l-/) =^ /(?/) = TT^- 

Therefore, one can expect a branching of TW profiles for small |A| > from (ll.40p at 
A = 0. Then, if this is not a subcritical pitchfork A-bifurcation, we expect existence f{y), 
at least, for all small A > (thought, not extensible to A = +1, as mentioned above). 

8. QuASILINEAR dispersion KPP-(11,1) PROBLEM 

Consider the quasilinear ODE in (ll.44p . with n > 0. First of all, we justify that the TW 
profiles are assumed to have finite interfaces at some finite point y = yo > 0; see extra 
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0.8 - 
0.6 
0.4 ■ 
0.2 - 









V ■ 








— X=O.S 










. 


m\ \ \ \ \ m 




f 1' V V y f f \i 



-30 -20 -10 10 20 30 40 50 

y 



(a) l^rm . A = 0.5,0.8,1 

f(y) KPP-(10,6):)u'f'°' = f''°' + f(1 -f),TW profile for). = 1 
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(c) (f03| . A = -1 



f(y)i 



KPP-(10,4): f'^'=- f''°'+f{1-f), TW profile for^ = 1 





(b) ((02l) . A = 1 

f(y) KPP-(10,8): f'"' = - f''°' + f(1-f), TW profiles for). = 1,2,3 




(d) (1131, A = 1,2,3 



Figure 10. TW profiles for hyperbolic equations (fr3T|) - pr3i|) . 



details for such a functional setting in (the singular conditions as ?/ — > — oo remain 
principally the same as in f ll.lOp ). It is clear that, as y ^ , the source term /(I — /) 
can be neglected, so that the asymptotics of f{y) near this finite interface is described by 
both leading higher-order terms, so one needs to study the following asymptotic ODE: 

(8.1) _A/' = -(|/r/)("^ =^ A/=(|/r/)^'°^ (A>0). 

Similar to the approach in [10] for quasilinear parabolic equations, we claim that this 
asymptotic behaviour is given by 

(8.2) f{y) = {yo - yVvi.s), where s = \n{yo - y) ^ -oo as y ^ y^ , 7 = ^> 
where the oscillatory component (/j(s) satisfies another complicated tenth-order ODE: 



(8.3) 



P-^^W^Y^ip] = \ip in M. 
15 




-300 -250 -200 -150 -100 -50 50 -2000 -1500 -1000 -500 




(c) (fOHll . A = 10 (d) (fOTl) . A = 10 

Figure 11. Oscillatory one-sided "periodic" TW patterns for elliptic equations 

(HMD-dOZD. 

Here, the linear polynomial operator Pio[0] belongs to the family {Pfc[0], A; > 0} of oper- 
ators that are constructed by the iteration 

(8.4) Pk+M = {PkW + (7 - for A; = 0, 1, ... , FqM = <P- 
In particular, this yields: 

(8.5) FsM = r + 3(7 - 1)0" + (37' - 67 + 2)0' + 7(7 - 1)(7 - 2)0; 



P4[0] =0(4) + 2(27 - 3)0"' + (67^ - I87 + 11)0" + 2(27^ - 97^ 
+ II7 - 3)0' + 7(7 - 1)(7 - 2)(7 - 3)0; 
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KPP-(10,11):-)." t'"'=-t"°' + f(1-f), TW profile for). = 0.2 
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KPP-(10,11 ):-)." f'"'= -f'"" + f(1-f),TW profile for). = -1 and -1.01 
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(a) arm , a = o (b) doli), a = -1,-1.1 

Figure 12. Oscillatory one-sided "periodic" TW patterns for elliptic equations 

([LlSD-dOTD. 
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KPP-(10,11): -A," f'"'= - f*^"' + f(1-f), TW profile forX = -1 
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(8.7) 



Figure 13. A TW profile for the KPP-(11,12), ffTaoD . for A = 1 and 0.8. 

P,[(l)] =0(5) + 5(7 - 2)(^(^) + 5 (272 - 87 + 7)0'" 

+ 5(7 - 2){2-f^ - 87 + 5)0" + (57^ - 407^ + 1057^ 

- IOO7 + 24)0' + 7(7 - 1)(7 - 2) (7 - 3) (7 - 4)0, etc. 



The operator Pio in (18. 3p is too ambiguous to present it here. 

The next main point now is that the ODE (18. 3p for admits a periodic solution (p^,{. 
which, together with its stability set as s — )■ —00 (including the obvious translations 
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(c) (fOSj) . A = -1 
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(d) (OjlL A = -1 
Figure 14. TW profiles for equations p^l0]) ~ pri3]) . 
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Figure 15. A TW profile F{y) satisfying (18. 8 p for A = 1 and n = 1. 

s), represents the actual asymptotic bundle of all admissible solutions of the form (18. 2p 
near the interface. 

Examples of such periodic solutions V5*(s) for operators Pq in (18. 3 p can be found in [HI 
p. 192]; see also similar examples for P5 in The higher-order case of Pg (the parabolic 
one) is given in |TT1 p. 143], etc. The Pio case is no much different, though, since (/?*(s) 
can be more unstable as s — )■ +00 (its unstability manifold s — )■ —00, i.e., approaching 
the interface, becomes even more dimensional), so numerics may get more difficult. 

Examples of full global solutions of the ODE in (18. ip are presented in a number of 
figures below, where, for convenience and by obvious reasons, we represent the function 

(8.8) F{y) = \f{yTf{y) =^ F(ii) = ^ + |Fr^F(l - 

Namely, Figure [15] shows the TW profile for A = 1 and n = 1. The next Figure [16] shows 
TW profiles F{y) again for n = 1 for A = 0.5 and A = 1.196. The last value turns out to 
be close to the maximal value Amax and we get the estimate 

(8.9) n = l: 1.196 < Amax(l) < 1-197. 

Recall that, for the semilinear case n = 0, Amax(O) is slightly larger; see (12. 2p . 

TW profiles for some negative velocities A and n = 1 are presented in Figure [17] For 
the sake of comparison, we also indicate therein the stationary profile, with A = 0, for 
n = 0. 

In Figure [T8| the TW profiles F{y) correspond to a larger n = 2 and A = 0, 1,2. The 
last value is not that far from the maximal value: our computations show that 

(8.10) n = 2: 2.25 < Amax(2) < 2.26. 
For n = 4, Figure [19] shows the profiles F for A = 0, 0.3, and 0.5. 
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F(y) KPP-(11,1)n: -Xf = -(|f|"f)*"'+f(1-f), TW F(y) forX=0,1,2; n=2 
1.2 




Figure 18. TW profiles F{y) satisfying (EHD for A = 0, 1, 2, and n = 2. 

F(y) KPP-(11,1)n: -Xr = -(|f|"f)<"'+f(1-f), TW F(y) forX=0,0.3,0.5; n=4 
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Figure 19. TW profiles F{y) satisfying (!8l8|) for A = (stationary), 0.3, 
0.5, and n = 4. 

9. A logt-SHIFT IN THE DISPERSION KPP-(11,3) PROBLEM: A CENTRE SUBSPACE 

PATTERN 

Thus, we begin with a formal analysis of a kind of a "centre subspace behaviour" , which 
generates a necessary logt-shift of the wave front. We restrict to a semilinear PDE. For 
similar applications to quasilinear (parabolic) ones, see |I0] . 

21 



As a typical PDE, bearing in mind fll.lSp . consider a semilinear KPP-type problem for 
a PDE, 



(9.1) uttt = A.u + u{l - u) in 



^+5 



where A is a proper homogeneous isotropic and translational invariant linear differential 
operator satisfying some extra conditions specified below. For instance, we can fix the 
dispersion operator 

(9.2) A = -Dl' (cf. dnZD). 
We assume that the corresponding ODE problem 

(9.3) -Xlr = Af + f{l-f), 

with the conditions (ll.lOp admits a unique solution /. 

Attaching the solution u{x,t) to the front moving and setting, as usual, Xf{t) = Aot — 
g{t), the PDE reads 

(9.4) u{x,t) = v{y,t), y = X - Xot + g{t). 
Then v satisfies the following perturbed equation: 

vttt + ^vtty{-Xo + g')' + Svtyy{-Xo + g'f + ^vtyg" 

= AV + V{1 -V)- Vyyy{-X0 + C/')^ " 5Vyy{-Xo + g^O" ' Vyg'"- 



Thus, again, as usual, we assume that g'{t) — )■ as t — ?■ +oo sufficiently fast, i.e., at 
least algebraically, so that 

(9.6) \g"{t)\'^\g'{t)\ and |/'(t)| < |/(t)|, etc. for t > 1. 
We next "linearize" (19. 5 p by setting 

(9.7) viy,t) = fiy) + wiy,t). 

Then, using (19. 3p . yields the following perturbed equation, where, according to (19. 6p . we 
keep the leading terms only: 

Wttt-^XoWtty + ^Xlwtyy + 3wtyg" = Bw - 3X1 f'g" - f'd'" - 

(9.8) „ 
where Bw = Aw + (1 — 2f)w + X^Wyyy. 

Note that here we face an essentially non-autonomously perturbed flow, so we cannot use 
advanced semigroup theory; cf. [15]. Instead, we will apply formal asymptotic expansion 
techniques. 

Now, assuming that, in this g{t)-movmg frame, there exists the convergence as in (II. 6p . 
so that w{t) — )■ as t — )■ +oo, one can see that, under the hypothesis (19. 6p . the leading 
non-autonomous perturbations in (19. 8 p are those of order 0{g"{t)), since the rest of the 
terms are negligible as t — > +oo. Therefore, one needs to balance these major terms, but 
then the actual behaviour of g{t) for t ^ 1 (and, hence, proper logt-shifts of the front) 
will depend on the next matching. 
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Thus, under the hypothesis (19. 6p . the only possible way to balance all the terms therein 
(including the quadratic one —w"^) for t ^ 1 is to assume the following expansion: 

(9.9) w{y,t) = g'{t)ij{y) + e{t)ip{y)... as t +oo, where \e{t)\ <^ \g'{t)\ 
is still an unknown coefficient. Substituting (19. 9p into (19. 8 p yields 

+ 5"V + 3(/y + 5V)(-Ao + ^?') 

+ 3{g"r + e'^"){Xl - 2\og' + {g'f) + 3{g"^' + e'^')g" + ... 

(9.10) = g'Bi) + eB^-{g'fi)'^ -2eg'^^lj 

+ - (/'" + g'r + e^"'){-\l + 3\lg' - SXoigT + (gT) 

- 3{f" + g'r + e^'K-Xo + 9')g" - if + q'^' + e^')g" + ... , 

where we have omitted some obviously negligible terms. 

Using (19. 6p and (19. 9 p in balancing first the leading terms of the order 0{g'{t)) yields 
the inhomogeneous equation for ip: 

(9.11) 0{g'{t)) ( = 0(i), see below) : B^^ - SA^/'" = 0. 
Then balancing the rest of the terms in (I9.10p requires 

(9.12) g"{t) ~ -{g'{t)f ~ ^(t), i.e., g{t) = klogt, g'{t) = f , g"{t) = -A, e{t) = ^. 
Then, we obtain the second inhomogeneous singular Sturm-Liouville problem for ip: 

(9.13) O(^) : Bip = k{3Xof" - SAg^ - /') + + SAg^'" - 3\of"'). 

Thus, the first simple asymptotic ODE in (19.120 gives the log t-dependence as in (II. 7p . 
Finally, we arrive at the following system for {ifj,ip}: 

B^|J = 3Xlf'\ 

Bip = k{3Xof" - 3Xlr - n + + - SAoD- 

Solving this system, with typical boundary conditions as in (ll.lOp . allows then continue 
the expansion of the solutions of (19.80 close to an "affine centre subspace" of B governed 
by the obvious (by translation) spectral pair 

(9.15) Ao = and V^o(2/) = /'(?/)• 
The asymptotic expansion for t ^ 1 then takes the form 

(9.16) w{y,t) = \i,{y) + j,p,{y) + ..., 

which can be easily extended by introducing further terms, with similar inhomogeneous 
Sturm-Liouville problems for the expansion coefficients. 

Since B does not have a discrete spectrum, we cannot derive a simple algebraic equa- 
tion for k by demanding the standard orthogonality of the right-hand side in the second 
equation in ( ]9.14p to the adjoint eigenvector -i/^Q of B* in the L^-metric (in which the 
adjoint operator B* is obtained), like 

(9.17) k : {k{3X,f" - 3Xlr - n + k\r + 3Ag^'" - SAo/'"), V^o) = 0- 
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(9.14) 



Therefore, the system fl9.14p cannot itself determine the actual value of k therein. As 
we have mentioned, the latter requires a difficult matching analysis in Inner and Outer 
Regions, which, for all present KPP-problems, remains an open problem. 

One can see that the above elementary conclusion well corresponds to a "centre subspace 
analysis" of the non-autonomous PDE (19.81) . and then r = logt naturally becomes the 
corresponding "slow" time variable; see various examples in [12] of such slow motion along 
centre subspaces in nonlinear parabolic problems with global and blow-up solutions. In 
the latter case, the slow time variable is r = — ln(T — t) — t- +oo as t approaches finite 
blow-up time T~ . 

For the semilinear higher-order reaction-absorption equations such as 

(9.18) Ut = —u^xxx — \u\^~^u in R x R_|_, with p > 1, 

existence of log t-perturbed global asymptotics was established in [S]. For finite-time 
extinction, with —1 < p < 1 in (I9.18p . this was done in [6]. For the corresponding 
blow-up problem with the combustion source 

Ut = —Uxxxx + |M|^~"^tt in M X ]R_|_, with p > 1, 

log(T — t)-dependent blow-up singularities were constructed in [4J. We must admit that 
any justification of such log t- front corrections (and finding corresponding classes of initial 
data) in any of KPP-(A;, /) problems with fc, / > 2 is much more difficult and remains an 
open problem. 

10. Final conclusions 

1. In all of KPP-(A;, /) problems, there exist TW solutions for various values of the 
speeds A. Most of them satisfy singular boundary conditions at ±oo. However, for some 
types of PDEs involved, these can be oscillatory and/or periodic either as y — -|-oo, or 
?/ — )■ — oo, and hence require special setting. 

2. In all of the higher-order KPP-(fc, I) problems, we did not observe the classic KPP- 
2 phenomenon for (II. ip . (II. 3p [T3] of existence of the minimal speed Aq = 2, such that 
TWs exist for A > 2 only. It seems that this phenomenon is directly connected with 
the Maximum Principle (and other features of order-preserving fiows), and becomes non- 
generic and non-existent when this fails. 

3. Moreover, in several KPP-(A;, /) problems, on the contrary, we observed existence 
of a maximal speed Amax, so that for slightly A > Amax, TW profiles do not exist. In 
particular, this always happens for parabolic problems, [9]. 

4. In all the KPP-(fc, I) problems, there exists a formal justification of existence of 
logt-drift of the propagating fronts in PDE setting along an "(affine) centre subspace" 
of semilinear rescaled operators involved. Then r = logt naturally appears as the corre- 
sponding "slow" time variable. 

5. We must admit that the important PDE problem on the actual structure of the 
omega-limit set uj{H) of the rescaled (properly shifter) solution u{x,t) of the various 
KPP-problems with the Heaviside data H{—x) remains open. In particular, it is not still 
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known whether uj{H) = {/(■; Aq)} for some Aq € A, i.e., whether uj{H) consists of a single 
TW profile. We believe that this problem deserves further study by analytical and PDE 
numerical methods, but expect it to be very difficult. 
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